^■■■■IHHHn 


A  NONLINEAR  CLUTCH  SERVOMECHANISM 
FOR  UNDERWATER  CONTROL  SURFACE 
ACTUATION 

by 

ROBERT  J.  CHENEY,  JR. 


DUDLEY  KNOX  LIBRAHY 
NAVAL  POSTGRADUATE  SCHOOL 
MONTFREY  CA   93943-5101 


A  NONLINEAR  CLUTCH  SERVOMECHANISM 
FOR  UNDERWATER  CONTROL  SURFACE  ACTUATION 

by 

ROBERT  J.    CHENEY  JR. 
/j 


B.S.     U.S.    Coast  Guard  Academy 
I960 


SUBMITTED  IN  PARTIAL  FULFILLMENT  OF  THE 
REQUIREMENTS  FOR  THE 
PROFESSIONAL  DEGREE  OF 
NAVAL  ENGINEER 
AND  THE  DEGREE  OF 
MASTER  OF  SCIENCE  IN  ELECTRICAL  ENGINEERING 

at  the 


MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 

August  1966 


A  NONLINEAR  CLUTCH  SERVOMECHANISM  FOR  UNDERWATER 
CONTROL  SURFACE  ACTUATION 

by 
ROBERT  J.  CHENEY,  JR. 

Submitted  to  the  Department  of  Naval  Architecture  and  Marine 
Engineering  and  the  Department  of  Electrical  Engineering  on  August  22, 
1966  in  partial  fulfillment  of  the   requirements  for  the  professional 
degree  of  Naval  Engineer   and     the     degree  of  Master  of  Science  in 
Electrical  Engineering. 

ABSTRACT 

A  nonlinear  clutch  type  servomechanism  is  proposed  to  position 
the  control  surfaces  of  an  underwater  vehicle  in  accordance  with  specified 
performance  criteria. 

Two  system  configurations  are  discussed  and  compared: 

(1)  A  two-state  system  employing  pulse  width  angular  velocity 
modulation  on  counter  rotating  clutches. 

(2)  A  three-state  system  employing  counter  rotating  clutches  and 
a  brake  with  error  sampling. 

Steady  state  load  vibration  will  be  present  in  each  case.     The 
degrading  effects  of  these  vibrations  are  developed  and  discussed. 

Clutch  power  dissipation  phenomena  are  discussed. 

Actuator  specifications  are  examined  in  light  of  overall  system 
dynamics . 

The  present  state  of  clutch  technology  is  examined  to  fulfill  system 
requirements  . 
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CHAPTER  I 

INTRODUCTION  AND  SUMMARY 

1 .  0    PROBLEM  STATEMENT 

Switching  amplifiers,    in  general,    offer  higher  efficiency  and 
higher  power  per  unit  weight  and  volume  than  do  continuous  amplifiers. 
This   research  investigated  the  application  of  clutch  type  actuators 
operated  in  a  switching  mode  to  the  feedback  position  control  of  the 
control  surfaces   of  an  underwater  vehicle.      The  vehicle  was  assumed 
to  be  operating  at  maximum  speed. 

1.  1     SYSTEM  CONFIGURATION 

A  self-propelled  underwater  vehicle  when  properly  designed  with 
neutral  stability  in  yaw-rate  will  exhibit  fastest  response  in  that  mode 
of  motion.      A  control  surface  position  servo  satisfying  the   requirements 
of  yaw-rate  tracking  will  be  adequate  for  roll  and  depth  control. 

A  linear  small  signal  model  for  the  yaw-rate  loop  is  indicated 
in  Fig.    1.  1(a). 

The  individual  transfer  functions  are 

^o  .  A(s+a) 


5 
o 

6 
o 

s(s+b)(s  +  c) 
2 

CO 

n 

6. 
i 

2   j     -?r              x      2 

S         +    Z>L  CO     S   +03 

where  K..K-,  are  preassigned  gain  factors 

A,  a,  b,  c      are  speed  dependent  vehicle  parameters 


co  =  natural  frequency  of  servo  actuator  system 

n  ' 

-1- 


-2 


*-k> 


8i 

Rudder 
Servo 

% 

Vehicle 
Dynamics 

*o 

Rate 
Gyro 

Fig.  1.1  (a)     Yaw-Rate  Loop 
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t,      =     damping  ratio 


4j  .    =     command  yaw  rate 


6.   =   command  rudder  deflection 
l 

v|j    =    vehicle  yaw  rate 

6     =  actual  rudder  deflection 
o 

A  linear  model  for  the  control  surface  position  control  system  is 
indicated  in  Fig.    1.  1(b).      The  control  surface  system  frequency  response 
has  been  specified  as  flat  to  7   cps  with     -3dB    point  at  27  cps.  The 

torque  amplifier    K-     must  deliver  up  to   1000  in     -lb     torque  at  a 
maximum  speed  of  100  RPM. 

1.2    ACTUATOR  CONFIGURATION 

The  design  of  a  linear  control  system  for  control  surface  positioning 
is  straightforward  and  can  be  realized  with  an  electric  or  hydraulic 
actuator.      These  devices  have  inherent  disadvantages,    notably  energy 
conversion  in  the  former  and  shelf  life  in  the  latter. 

Clutch  type  actuators   offer  superior  performance  in  many  respects 
although  analysis  is  difficult  due  to  nonlinear  describing  equations. 
Counter  rotating  clutch  actuators  can  be  operated  in  a  push-pull  fashion 
or  in  a  switching  mode.      The  former  is  undesirable  from  efficiency 
and  heat  sinking  considerations.     Switching  operation  is   considered  in 
this  research. 

Two  basic  configurations  are  considered:     A  two-state  pulse  width 
modulated  counter-rotating  clutch  arrangement  and  a  three-state 
system  utilizing  two  counter-rotating  clutches  and  a  brake  with  state 
determined  by  sampled  errors. 
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The  two-state  system  is  functionally  diagrammed  in  Fig.    1.2(a). 
The  pulse  width  modulator  output  is  designed  to  operate  at  a 

constant  recurrence   rate  with  fixed  output  magnitude  positive  or  negative 

(3) 

with  half  cycle  duration  determined  by  input  level.  The  pulse  width 

modulator  output  drives  the  actuation  mechanism  of  the  clutch. 

High  inertia  clutch  input  members  are  driven  at  constant  speed 
from  the  vehicle  prime  mover.     Assume   that    input  and  output  members 
become  instantly  locked  on  signal  from  the  pulse  width  modulator. 
The  clutch  output  member  velocity  will  be  a  square  wave.     The  spring 
Kf    transforms  this  velocity  square  wave  into  a  triangular  torque  wave 
whose  short-time  mean  level  is  transmitted  to  the  load.     The  load  and 
filter  inertia  attenuate  torque  components  at  switching  frequency 
and  higher  harmonics. 

Compensation  and  feedback  networks  insure  reset  action  and 
stabilize  the  system.     Reduction  gearing  is  necessary  to  reduce  clutch 
size. 

In  the  steady  state  the  clutch  output  velocity  will  be  a  symmetrical 
square  wave.     Some  wobble  will  be  inherent  in  the  control  surface  load 
at  switching  frequency  and  its  harmonics. 

The  three  state  system  operation  is  self  explanatory  in  Fig.  1.2(b). 
System  operation  is  quite  similar  to  that  of  the  two-state  configuration. 
The  notable  exceptions  are  the  more  complicated  decision  box  and  lack 
of  steady  state  output  oscillations  at  the  expense  of  steady  state  errors,  e 
Limit  cycle  oscillations  can  occur  at  sampling  frequency  if  the  forward 
loop  gain  is  high  enough. 
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Detailed  design  of  the  clutch  actuator  system  is  the  subject  of 
a  separate  thesis  and  will  not  be  dealt  with  here. 

The  advantages  of  clutch  type  actuators  derive  from  the  high 
torque  squared-to-inertia  ratios  and  low  control  power  requirements. 
Shelf  life  is  no  problem  with  clutch  servos. 

The  disadvantages  of  both  configurations  include  heat  dissipation 
at  contact  surfaces  and  output  wobble  with  the  possibility  of  exciting 
structural  resonance. 

Selection  of  the  optimum   switching  of  sampling  frequency  will 
involve  tradeoffs  between  output  wobble  amplitude,    which  will  decrease 
with  switching  frequency,    and  power  dissipation  in  the  clutches,    which 
will  increase  with  switching  frequency.     The  fundamental  switching 
frequency  must  lie  above  load  resonance  and  away  from  structural 
resonances . 

1 .  3    CONTROL  SURFACE 

The  adverse  effects  of  control  surface  wobble  were  investigated 
and  acceptability  criteria  were  assigned. 

Control  surface  wobble  manifests  itself  as  a  perturbation  on 
the  control  surface  angle  of  attack.      The  hydrodynamics   of  the  situation 
were  examined.     Since  mean  incidence  remains  unchanged  the  dynamic 
lift  of  the  control  surface  is  not  affected.     Viscous  drag  effects  were 
second  order  in  magnitude,    in  fact,    with  the  proper  phasing  and  correct 
relationship  between  frequency,    vehicle  velocity  and  control  surface 
chord  length,    propulsion  will  result.      Cavitation  inception  was   considered 
unlikely  since  the  disturbance  velocities  were  insignificant  compared 
with  quasi-steady  levels. 
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Sound  generation  by  the  oscillating  control  surface  proved 
to  be  a  significant  result.      The  sonic  intensity  level  can  be 
determined  from  the  equation 


,dL.2         2 

(  *  *    C°5     °  „    /     2 

l(r)    =      2 5 T~        watts/m 

16  tt    p     c        r 
ro    o 

where  L    -      lift  force  exerted  by  hydrofoil 

0    ~      angle  between  direction  of  observation 
and  lift  force 

p      -      density  of  undisturbed  fluid 
ro  ' 

c       =      speed  of  sound  in  undisturbed  fluid 
o 

r     =     distance  from  source  to  observer 

This  relationship  is  developed  in  Chapter  Two.     Thus,    knowing 
lift  as  a  function  of  angle  of  attack  and  angle  of  attack  as  a  function 
of  time,    the  sonic  intensity  can  be  obtained.      The  sound  pressure   level  (SPL) 
can  be  obtained  from  intensity.     A  convenient  measure  of  sound  is  the 

source  level  (far  field  radiation  extrapolated  to     r  =   1M)  expressed  in 

-  5  2 

decibel  ratio  to  0.002  microbar  (2  x   10         Newton/Meter    ). 

For  a  given  maximum  control  surface  rate  specification  (100  RPM 
in  this  case)  the  source  SPL  can  be  determined  as  a  function  of  switching 
frequency  in  the  two-state  configuration.     A  limit  cycle  amplitude 
determination  must  be  made  for  the  three-state  configuration. 

The  two-state  source  SPL  was  observed  to  decrease  20  dB  per 
decade  of  switching  frequency. 
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1.4    CLUTCH  POWER  DISSIPATION 

Referring  to  Fig .    1.2(a)  if    Kf    were  zero,  the  amount  of  energy 

necessary  to  accelerate  the  clutch  output  member  from  full  clockwise 

angular  velocity  is 

E  ,     =     J    f2      joules 
d  c  J 

J         =     clutch  output  inertia 

c  r 

fi        =     angular  velocity  of  input  members 

Since  this  operation  is  performed  twice  in  a  switching  period, 

P,    =     2J  Q     f  watts 

d  c  sw 

f        =     switching  freq  (cps) 

The  dissipation  rating  of  the  clutch  must  lie  above  this  level. 
The  situation  can  be  alleviated  somewhat  by  making    K,    as  stiff  as 
possible  without  exceeding  clutch  pullout  torque  in  spring  windup.      An 
analog  simulation  bears  this  out  (see  Fig.    4.  1(b)  ).     This  result  is 
consistent  with  two-state  system  design  objectives  in  that  clutch 
pullout  torque  must  be  higher  than  the  greatest  demanded. 

In  the  analog  simulation  the  filter  inertia     J,    (Fig.    1.2(a))  was 

considered  infinite.     In  the  actual  system     Jr    >    5J       for  efficient 

3  f  c 

filtering  so  this  assumption  was  not  a  bad  one.     The  simulation  is 
detailed  in  Chapter  Four. 

Since  power  dissipation  is  a  linear  function  of  switching 
frequency  and  a  quadratic  function  of  input  member  speed  the  gear  ratio 
should  be  selected  to  minimize  input  member  speed.      Unfortunately 
load  torque  requirements  determine  the  gear  ratio  since  the  clutch 
dimensions  must  be  kept  small.      Pullout  torque  is  determined  by  clutch  size, 
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If  the  control  surface  can  be  pivoted  nearer  to  its  center  of 
pressure  the  output  torque  required  (1000  in -lbs)  will  be  reduced 
and  the  necessary  reduction  ratio  can  be  reduced  permitting  lower 
clutch  input  member  speed  while  maintaining  maximum  control  surface 
position  rate  at  the  specified  100  RPM. 

1.  5    OVERALL  SYSTEM  DYNAMICS 

Since  the  control  surface  servo  is  embedded  within  the  overall 
vehicle  dynamic  control  loop,    a  study  was  carried  out  to  determine 
the  effect  of  actuator  bandwidth  reduction  on  the  dynamics  of  the 
total  system. 

The  vehicle  exhibits  fastest  behavior  in  the  yaw  rate  loop  (Fig.    1.  1(a)), 
The  control  surface  servo  system  was  modelled  as  a  second  order 
system  with  damping  ratio  equal  to    y  . 
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It  can  be  shown  that  damping  ratio  =0.5    minimizes  integral  square 
error  for  a  linear  position  servo  for  a  step  input.      This  behavior  can 
be  achieved  by  proper  selection  of  forward  loop  gain  and  output  rate 
feedback. 

Bode  plots   of  the  open  yaw-rate  loop  were  constructed  for  the 
specified  servo  and  three  slower  servos.     These  data  were  transferred 
to  a  Nichols  chart  for  comparison  (Fig.    1.5). 
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Fig.  1.5     Nichols  Charr  of  Yaw-Rate  Loop   Frequency  Response 
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It  was  concluded  that  the  response   speed  of  the  control  surface 
actuator   could  be  halved  without  degrading  the  closed  yaw-rate 
loop  performance. 

1.6    CLUTCH  SPECIFICATIONS 

The  state  of  clutch  technology  was  scanned  and  unfortunately 
present  production  line  clutches  do  not  fulfill  system  requirements. 
Switching  rates  in  the  range  75  to   150  cps  are  necessary  for  the  two 
state  system.      This  gives   rise  to  an  input-output  rise  time  of  one  or 
two  milliseconds.      This  demand  has  been  met  by  very  small  instrument 
magnetic  particle  clutches.     Figure   1.6  compares  characteristics  of  three 
production  devices  with  proper  torque  and  speed  characteristics. 

Maximum     Maximum      Maximum       Rise  *\.  Dissipation 

Clutch        Torque  Speed  Dimension    Time  _  .     Rating 

Type  (in-lb)  (RPM)  (in)  (±n  sec)    (in-lb  sec   xlO    )    (watts) 


Solenoid- 
Operated 
Mech. 
Friction 

Magnetic 
Particle 


60 


35 


7500 


3400 


3.7 


3.  1 


30 


30 


6.93 


2.37 


136 


75 


Solenoid- 
Actuated 
Vacuum-        60 
Operated 
Mechanical 


1800 


9.3 


127 


Fig.    1.6    Comparative  Clutch  Characteristics 
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In  addition,    the  magnetic  particle  clutch  exhibits   0.6  in-lb  drag 
torque  at  zero  excitation  and  the  vacuum -operated  clutch  requires   a 
small  vacuum  pump. 

It  may  be  concluded  that  the  two-state  actuator  configuration  is 
not  technically  feasible  at  present.     Since  the  three-state  system  has 
not  been  fully  investigated,    a  working  model  may  be  constructed  using 
the  present  clutch  technology.     High  speed  magnetic  tape  drives  have 
indicated  a  need  for  fast  mechanical  clutches  and  current  research  may 
evolve  a  device  capable  of  meeting  the  requirements  of  the  two-state 
actuator  system. 


CHAPTER  II 
SOUND  GENERATION 

The  intensity  of  sound  resulting  from  a  lifting  body  in  fluid  flow 

(4) 
has  been  investigated.  Rayleigh  has   shown  that  any  frequencies   in 

a  fluid  flow  are  equal  to  those  of  the  sound  generated. 

Lighthill  has   considered  the  equations   of  motion  of  fluids  and 

compared  them  with  the  equations   of  sound  propagation  in  a  medium 

at  rest.      The  exact  equations  for  the  former  may  be  written  in  the  form 

A   -   c2   v2    =  r4-    T..  (2.D 

where  T  .  •      =        p  v.v.  +  p.  .   -   C       6.. 

ij  i    J  ij  o       lj 

p       =       density 

p..       =       compressive  stress   tensor 
ij 

C         =       sonic  velocity  in  fluid  at  rest 
o  ; 

v.       =      velocity  componentin   direction  x.      (i  =   1,  2,  3) 

Physically,    this  means  that  sound  is  generated  by  a  fluid  flow 
exactly  as  in  a  uniform  medium  at  rest  which  is   acted  upon  by  externally 
applied  fluctuating  stresses.     If  we  regard  the  flow  as  known  it  is  possible 
to  evaluate  the  sound  generated. 

Equation  (2.  1)  is  the  non-homogeneous  wave  equation.      The  most 
general  solution  (Stratton  1941)  is 


P(x>    "  PQ 


3r_   |£    dS(D 

n    at  w  ' 
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In  this  equation    — — *  ,     —  -  ■  ,    p   ,    ~ir~-     are  taken  at  retarded  times 

oy.dy.  3n  9t 

1      J 

t  -    —     ,    where     r  =    I  x  -   y~  I      and     rT  is   the  outward  normal  from  the 

c  '  '    ' 

o 

fluid.       The  first  integral  is  taken  over   the  total  volume     V     external 
to  the  solid  boundaries  and  the  second  integral  is  taken  over  the 
surface    S     of  the  solid  boundaries. 

Curie  has  shown  that  the  solution  may  be  simplified  under  the 
assumption  that: 

(1)  |  x*"|     »    X.        ("far  field"  approximation) 

(2)  the  solid  boundaries   are  rigid,    i.e.,    non -dilatational 

(3)  the  flow  is  incompressible  (Mach  No.    «   1) 
Under  these  assumptions  the  solution  reduces  to 

3      z       at     ili; 


P   -  P0    =     9     = 

4  ttc         x 
o 


where  F.(t)    =        /       P.(y  ,  t)  d  S(y  ),      the  total  resultant  force 

S 

exerted  on  the  fluid  by  the  solid  boundaries. 

In  acoustic  media,    the  sound  pressure  field  is  proportional  to 

2 

the  density  fluctuation  field  by  the  constant     c        thus 

1  '  o 

a  i        xi      a 


p    =     4Tc~    - z   at  Fi(t) 


o       X 

This  is   the  radiation  field  of  a  dipole  source  with  axis  parallel 
to     F(t)     the  net  resultant  force  exerted  on  the  fluid  by  the  boundaries 
Rewriting  in  more  convenient  variables; 
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A  1  cos    6      9     tp  ...  .,    _. 

P    =      ~a ~zz    F.(t)  (2.2) 

o 

where     0    is  the  angular  distance  between  the  radius  vector     r     and  the 

axis  of    F.(t). 
l 

Since  intensity  is  defined  as  the  acoustic  power  transmitted  through 
a  unit  area 


I 


2 
P 


P    c 
ro    o 


2  2 

p        =      mean  square  acoustic  pressure  (Newtons/M    ) 

p    c        =      characteristic  impedance  of  acoustic  medium 
°   °  (MKS  RAYLS) 

Generally  it  is  more  physically  meaningful  to  consider  a  decibel 
ratio  referred  to  a  datum  level  for  intensity  and  pressure  measurements. 
A  customary  quantitative  measurement  in  underwater  acoustics  is 
the  source  level  sound  pressure  level.     This  is  defined  as 

~~ 2 

SPL    =     10  log1Q      -2j-      dB      ' 

Po 

~T 

where    p       is  far-field  acoustic  mean  square  pressure  extrapolated 

to  a  radial  distance  of  one  meter  from  the  source,    the  reference  pressure 
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being  0.0002  microbar  (2  x   10      N/M    )     for  underwater  noise  measurements 

If  the  force  exerted  by  the  body  on  the  fluid  can  be  determined  then 

the  SPL  can  be  calculated. 


CHAPTER  III 
HYDRODYNAMICS  AND  CONTROL  SURFACE  OSCILLATION 

Principles   of  incompressible,    non-cavitating  hydrodynamics 
can  be  utilized  to  determine  the  external  forces  acting  upon  the 
fluid  which  will  yield  the  sound  pressure  field.     Sound  generated  by 
the  oscillating  control  surface  only  is  considered.      This  sound  pressure 
will  be  additive  to  any  presently  inherent  in  vehicle  operation  exclusive 
of  the  control  surface  system. 

In  Equation  (2.2),      F.(t)     must  be  determined  to  find  the  SPL.      This 
total  force  exerted  by  the  body  on  the  fluid  is  due  primarily  to  dynamic 
lift  directed  normal  to  the  approach  flow.      A  secondary  force  is  due  to 
viscous  drag  which  is  oriented  parallel  to  the  approach  flow  and  generally 
equal  to  about  one  twentieth  of  the  lift  in  magnitude. 

For  an  oscillating  hydrofoil  the  oscillating  component  of  lift  is  a 
function  of  approach  velocity,    hydrofoil  geometry,    fluid  density,    and  the 
magnitude  and  frequency  of  the  oscillation.     For  a  steady  hydrofoil  the 
dynamic   lift  is  a  function  only  of  approach  velocity,    foil  geometry, 

fluid  density  and  angle  of  incidence  with  approach  flow.     A  frequently 

(5) 
used  dimensionless  parameter  is   the  lift  coefficient     C     . 

L 

CT      =       ^--y    =      Kq        (-17°<a<17°) 

L  p  U2 

for  a  symmetrical  thickness  form. 
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where  L  =  dynamic  lift 

p  =  fluid  density 

S  =  area  of  foil  in  the  span-chord  plane 

U  =  approach  flow  velocity 

a  =  radian  angle  of  attack 

K  =  proportionality  constant  dependent  on  aspect  ratio 

The  aspect  ratio  is  the  ratio  of  span  squared  to  foil  area    S.     For  two 
dimensional  flow  (infinite  aspect  ratio),     K  =  2  tt.     For  zero  aspect  ratio 
(chord  length  much  greater  than  span),      K  =   tt/2  .     It  is  known  that  at 
reduced  frequencies  greater  than  0.5  the  oscillating   component  of 
the  flow  becomes  two-dimensional. 

K     =      -jty      =     reduced  frequency 

to        =      radian  frequency  of  oscillation 

Q 

y        -      semi-chord  length 
U       =      approach  flow  velocity 

In  the  course  of  a  mathematical  investigation  of  wing  flutter 
phenomena,    Theodorsen        has   shown  that 

L  H(,2)  (k) 

OSC        =    C(k)    =      - 


Lstdy  H(2)  (k)  +  jH^2)  (k) 


for  a  lifting  surface  in  rotary  oscillation. 

where  L  =       complex  value  of  oscillatory  lift  component 

osc  ^  ' 

L    .  ,      =       quasi-steady  lift  corresponding  to 


maximum  inclination 


r<2)      - 


H  =       Hankel  function  of  second  kind  of  order     n 

n 


-1 


The  function    C(k)    is  plotted  in  Fig.    3.  1.     Its  behavior  is  quite 
similar  to  that  of  a  lag  network.      The  reduced  frequencies   of  the 
vehicle  control  surface  in  the  range  of  interest  were  found  to  be  greater 
than   1 .  0  so  that     C(k)    was  assumed  equal  to      [  0.  5  +  j  0.  0]     throughout 
the  analysis . 

The  foregoing  conclusions  have  been  predicated  on  the  assumption 
that  the  foil  is  oscillating  in  a  sinusoidal  fashion  with  time. 

Referring  to  Fig.    1.2(a),    assume  that  the  clutch  is  perfect  in  that 
the  output  member   locks   on  with  the  appropriate  input  member  at  the 
switching  instant.      The  angular  velocity  of  the  clutch  output  member  will 
be  a  square  wave  with  half-cycle  duration  determined  by  the  pulse  width 
modulator  input.      A  Fourier  decomposition  of  this   signal  yields: 


co 


jmto   t 


ZJ1I1UI 
F(m)  e 

m  /  0 


m  =  -  co 


-jmco   t 
2A(l-e  1    ) 

jmco 


F(m)    =< 


1 


(2  1   -    1)A      ;     m  =  0 


where 


Q 


to 


1 


=      clutch  output  member  angular  velocity 

=      switching  rate  (rad/sec) 

-      magnitude  of  input  member  speed 


t         =      duration  of  positive  half-cycle  of  PWM  output 

m       —      harmonic  index 

T       ~      period  of  PWM  output 
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In  the  case  of  steady  state  zero  position  error,     t  = 


The 


magnitude  of  the  fundamental  component  of  clutch  output  velocity  will 

4A 
be  a  maximum  of    . 


to 


Referring  again  to  Fig.    1.2(a)  relations  between  control  surface 
output  velocity  and  clutch  velocity  can  be  written: 


c 


a 


7  2T 

s  Z  s    J 


where 


L 
n 


load  impedance 
reduction  ratio 
s       =       complex  frequency 

For  switching  frequencies  above  load  resonance  (25  cps),    inertia 


terms   dominate  and    Z 


Recombining:      — —    =     n 
S2L 


1     + 


s2(Jf+    -^)n 
n 


Examining  the  relative  magnitudes  of  inertias: 

-3  2 

Jf     x      3x10        in  lb  sec     (3-5  J„  ,    see  Fig.    1.6) 

-2  2 

J       =      8x10        in  lb  sec 

Li 

since  a  reduction  ratio  of  at  least  15:1  will  be  required,    the  effects 
of  reflected  load  inertia  can  be  ignored  with  the  result  that 


L 


1 


a 


2     Jf 
n[l  +s     -jf] 


indicating  that  switching  frequency  harmonics  will  be  attenuated 


inversely  as  the  square  of  the  order.      Coupling  this   result  with  the    — 
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character  of  the  clutch  output  spectrum,    the  foil  motion  can  be 
considered  a  pure  sinusoid  at  fundamental  switching  frequency. 
The  source  SPL  can  now  be  calculated 


_dL   2 

v     dt    ' 


source 


1A      2         2 
lb   TT      c 

o 


where  L    =     L  sin  co,t    =     C(k)    KT    \9T    sinco,t1 

osc  1  L  L    L  1    J 


KT      =     constant  of  lift-incidence  proportionality 
K        =     TrpS   U2     (MKS  units) 
C(k)    =    Theodorsen's  function 


Since  the  sound  pressure  is  a  result  of  vortex  shedding,    a 
measure  of  spanwise  phase  correlation  is  necessary.     A  good  measure 
of  phase  correlation  is  obtained  by  a  strip  theory  approach.     The  foil 
was  arbitrarily  subdivided  into  four  spanwise  sections.     The  Theodorsen 
function  was  found  to  be  nearly  invariant  across  the  span  so  that  the 
sound  field  was  assumed  100  percent  phase  correlated.     The  tip  vortex 
baffle  increases  the  reliability  of  this  assumption. 

Converting  the  mean  square  pressure  to  decibel  ratio 

SPL  =     20  1og.„(^T    )    +    P.  dB 

source  b  1 0  v     L.  1 


=     20  Lot 


10 


2     Jf 
Ln(l-co1-^-). 


+    P        dB 


P.     contains     K       ,    C(k),    reference  pressure  (0.002     |xb),    (0.707)    , 
1  J-j 


16  tt       and     C 
o 
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Q 

Q 

If  we  further  assume  that    =     100  RPM     (specified  maximum 

n  v    r 

control  surface  speed),    then 

SPL  =     P,   -20  1ogin      o>   (1-u)2-     -~-    ) 

source  2  b 1 0         lv  1      Kf    ' 

This   relationship  is  plotted  in  Fig.    3.2  for  various   ratios  of 


Kf 
switching  frequency  to     / — = — 


/ 
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Fig.  3.2     Source   SPL  vs   Switching   Frequency 


CHAPTER  IV 
CLUTCH  POWER  DISSIPATION 

As  mentioned  in  Section  1.4  power  is  dissipated  in  the  interval 
between  the  switching  instant  and  the  point  at  which  clutch  input  and 
output  members  become  locked.      The  instantaneous  dissipation  is  the 
product  of  pullout  torque  and  slip  speed.     The  time  average  power 
dissipated  will  be  approximately  twice  the  product  of  clutch  output 
inertia,    input  member  speed  squared  and  switching  rate. 

In  order  to  minimize  power  dissipation,    the  absolute  value  of  slip 
velocity  must  be  minimized.      Assuming     J     ,    J,    and  switching  frequency 
fixed  (see  Fig.    l.Z(s)),    the  effect  on  slip  velocity  of  varying    Kf    was 
observed  on  an  analog  simulation.     Since     J-    will  be  three  to  five 
times  greater  than     J     ,    J,    was  assumed  infinite  for  purposes   of 
simulation.     The  system  simulated  and  computing  diagram  are  indicated 
in  Figure  4.  1(a). 

The  simulation  was   carried  out  on  the   Philbrick  Analog  Computer 
at  the  M.I.T.    Engineering  Projects  Laboratory.     Results  are  indicated 
in  Figure  4.  1(b).      The  input  member  velocity  is  the  square  wave. 
The  output  member  velocity  is  the  initially  ramp-like  square  wave 
which  tracks  the  input.      The  third  trace  is  the  amplified  difference 
between  the  two.     A  small  biasing  error  is  evident  in  seme  of  the 
computations.     This  is  attributed  to  amplifier  drift  since  it  was  necessary 
to  saturate  one  amplifier  in  the  loop. 

From  the  simulation,    it  may  be  concluded  that  the  spring    Kf 
should  be  made  as  stiff  as  possible  without  exceeding  clutch  pullout 
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Jc 


Kf 


1 


SYSTEM    SIMULATED 


TM 

-TM— 

COMPUTING  DIAGRAM 


Fig.  4.1  (a)     Analog  Simulation  of  Clutch  Dissipation  Phenomena 
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=    1.21 
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=    0.735 
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00 
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00 
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=    1.385 


Note  slip  as  spring  is  wound 
up  to  pullout  torque 
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00 


=    1.0 


sw 


00 


00 


=    1.54 


SW 


00 


u         =  switching  freq. 


SW 


Note   slip  as   spring  is  wound 
up  to  pullout  torque 


Fig.    4.  1(b)     Response  of  Second  Order   Clutch  Output  System 


-27- 
torque  in  spring  windup.     At  first  glance,    no  differences  may- 
be apparent  in  Fig.    4.  1(b).     A  closer  examination  reveals  that  response 

o>  co 

time  for    ~JL~      -   1.21     is  about  half  that  for      -2—     =     0.478. 

CO  CO 

sw  sw 

Input  member  speed  and  switching  rate  must  be  accounted  for 
in  the  final  selection  of    K,.. 
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